arXiv: 1508.01183v2 [math.GT] 31 Jan 2016 


LINKING NUMBER AND WRITHE IN RANDOM LINEAR 
EMBEDDINGS OF GRAPHS 

ERICA FLAPAN AND KENJI KOZAI 


Abstract. In order to model entanglements of polymers in a confined 
region, we consider the linking nnmbers and writhes of cycles in random 
linear embeddings of complete graphs in a cube. Our main results are 
that for a random linear embedding of Kn in a cube, the mean sum of 
squared linking numbers and the mean sum of squared writhes are of the 
order of 9{n{n\)). We obtain a similar result for the mean sum of squared 
linking numbers in linear embeddings of graphs on n vertices, such that 
for any pair of vertices, the probability that they are connected by an 
edge is p. We also obtain experimental results about the distribution of 
linking numbers for random linear embeddings of these graphs. Finally, 
we estimate the probability of specific linking configurations occurring 
in random linear embeddings of the graphs and it's, 3 , 1 . 


1. Introduction 

Long polymers become tangled up as a result of being tightly packed in a 
confined region. For example, 46 human chromosomes are packed together 
inside the nucleus of a cell whose diameter can be as little as 10“® times 
the length of a single chromosome. The entanglement that results affects 
the processes of replication and transcription of the DNA. For synthetic 
polymers, tangling is correlated with viscoelastic properties, and hence is 
important in the design and synthesis of new elastic materials. As these 
examples illustrate, understanding the tangling of polymers is useful for 
explaining and controlling molecular behaviour. However, since detailed 
visualizations of molecular entanglements are not yet technologically possi¬ 
ble, their study has been approached through mathematical modeling rather 
than experimental observation. 

Many authors have considered uniform random distributions of open and 
closed polygonal chains in a cube as a model for long molecular chains 
in a conhned region (see for example m, a, IS], 0, m, m, [is]). of 
particular note, Arsuaga et al |1] obtained a formula for the mean squared 
linking number of two uniform random n-gons in a cube, and showed that the 
probability of linking between a given simple closed curve in the cube and a 
uniform random n-gon grows at a rate of at least 1 — 0 • More recently, 
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Panagiotou et al m has shown that the mean squared linking number, the 
mean squared writhe, and the mean squared self-linking number of oriented 
uniform random open or closed chains with re vertices in a cube all grow at 
a rate of 0{n?). 

Many of the results that have been obtained about the linking of two 
random uniform polygons in a cube are restricted to pairs of polygons of 
the same length. However, there is no biological reason for all polymers in a 
given region to have the same length. While no theoretical results have been 
proven thus far about linking between uniform random re- and m-gons in a 
confined region, Arsuaga et al [T] has observed from numerical simulations 
that the linking probability of a random linear re-gon and rre-gon in a cube 

seems to be bounded below by 1 — O 


In order to obtain theoretical results about linking between random chains 
of different lengths as well as to measure entanglement in a more general 
way, we take a new approach. In particular, the above models use ordered 
sequences of re points chosen from a uniform random distribution of points 
in a cube to define one or two re-gons with linear edges. By contrast, we 
begin with an unordered set of re points chosen from a uniform random 
distribution of points in a cube. We then consider every possible pair of 
disjoint polygons obtained by adding line segments between some number 
of points in the set. By taking the sum of the squared linking numbers of 
all such pairs of polygons, we obtain a single number which represents the 
linking of all pairs of polygons with vertices in this set regardless of whether 
the polygons are of the same or distinct lengths. In addition, in order to 
measure the entanglement of individual polymers in a confined region, we 
consider the sum of the mean squared writhes over all polygons with vertices 
in our set. 

In addition to modeling entanglement of polymers in confined regions, 
our results can be seen in the context of linear embeddings of graphs in 
(that is, embeddings whose edges are realized by straight line segments). In 
particular, the set of polygons we are considering are the cycles in a linear 
embedding of the complete graph Kn in R^. Probably the most significant 
result in the study of such embeddings was the proof by Negami m that 
for every knot or link J, there is an integer R{J) such that every linear 
embedding of the complete graph Kr(j) in R^ contains J. In addition, 
several authors have obtained results characterizing what links can occur 
in linear embeddings of specific graphs. In particular, Hughes [7] and Huh 
and Jeon [ 8 ] gave combinatorial proofs that every linear embedding of Kq 
contains either one or three Hopf links and no other links. More recently, 
Nikkuni [12] obtained the same result with a topological proof. Naimi and 
Pavelescu [9] proved that every linear embedding of Kg contains a non-split 
link of three components, and showed in m that every linear embedding 
of A’ 3 ^ 3^1 contains either 1,2,3,4, or 5 non-trivial links. 
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Our main results concern the rate of growth of two measures of entangle¬ 
ment. Since rates of growth can be measured in several ways, for clarity we 
make the following definitions. 

Definition 1.1. Let /(n) be a function of the naturals. 

• /(n) is said to be of the order of 0{g{n)) if there exists a constant 
C > 0 such that for sufficiently large n, 

/(n) < Cg{n). 

• /(n) is said to be of the order of 9{g{n)), if there exist constants c, 
C > 0 such that for sufficiently large n, 

cg{n) < f{n) < Cg{n). 

Section 2 is devoted to the proofs of the following two theorems about 
entanglement of random linear embeddings of complete graphs inside a cube 
= [0,1]^. That is, embeddings of complete graphs whose vertices are 
given by a random uniform distribution of n points in the cube and whose 
edges are realized by straight line segments. 

Theorem 12.41 Let n > 6, and let he a random linear embedding of the 
complete graph on n vertices in the cube . Then the mean sum of squared 
linking numbers for is of the order of 6{n{n\)). 

Theorem 12.61 Let n > 3, and let he a random linear embedding of the 
complete graph on n vertices in the cube . Then the mean sum of squared 
writhe for Kn is of the order of 9{n{n\)). 

The ideas of the proofs of these results are as follows. The first theorem is 
proved using Lemma 12.31 which shows that the expected linking number of 
random cycles of length k and length I is on the order ol 9[kl). The complete 
graph Kn contains on the order of 9[{n — 1)!) links, both of whose cycles are 
of length Heuristically, the linking number from these cycles dominate 
the sum of squared linking numbers, resulting in a mean sum of squared 
linking numbers on the order of 9(n{n\)). To prove the second theorem, 
we observe that there are many more cycles of length n than of any other 
length, and the number of such cycles is on the order of 9{{n — 1)!). We 
then use the result of Panagiotou et al m that the mean squared writhe of 
an n-cycle is on the order of 9(n‘^) to obtain the desired result. 

In Section 3, we consider a set of n points chosen from a uniform random 
distribution of points in a cube, and then assign a probability that a given 
pair of vertices is joined with an edge. In this way, we obtain a subgraph of 
Kn with a given probability. In particular, we use the following definition, 
originally due to Gilbert [6]. 

Definition 1.2. Let n G N and p G (0,1). Then a {n,p)-graph is a graph 
on n vertices, such that for any pair of vertices, the probability that they 
are connected by an edge is p. 
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We obtain several results on knotting and linking in random linear em¬ 
beddings of (n,p)-graplis, including the following. 

Theorem 13.IL For any knot J, the probability that a random linear embed¬ 
ding of an {n,p)-graph contains a cycle isotopic to J goes to 1 as n ^ oo. 

Theorem 13.31 For any n > 6 and p G (0,1), the mean sum of squared 
linking numbers of a random linear embedding of an {n,p)-graph is of the 
order of 6{p^n{n\)). 

In Section 4, we apply the inequalities that we obtain in Section 2 to 
random linear embeddings of the graphs Kq and K 3 , 3 ,i in order to estimate 
the probability of specific linking configurations occurring. 

Finally, in Section 5, we sample random linear graph embeddings and 
make some observations about the distribution of linking numbers. 


2. Random linear embeddings of Kn 

Definition 2.1. A random linear embedding of a graph G, is an embedding 
of G in the unit cube G^ = [0,1]^ such that the vertices of G are embedded 
with a uniform distribution, and every edge {vi,Vj) of G is realized by a 
straight line segment between Vi and Vj. 

Arsuaga et al [T] prove the following lemma. 

Lemma 2.2 ([I]). Let Ii,l 2 ,li,l 2 denote edges in a random linear embedding 
of a graph with an orientation assigned to each edge. Let Cj denote the signed 
crossing of k and l[, and let E\eie 2 ] denote the expected value of eie 2 - 

(1) If the endpoints o/Zi, Z 2 , ^2 distinct, then E[eie 2 ] = 0. 

(2) If li = I 2 , the endpoints of l'^ and I 2 are distinct, and both are 
disjoint from Zi = I 2 , then E[ei€ 2 ] = 0. 

(3) Define variables as follows: 

• Let 2s denote the probability that Zi and l[ cross when Zi and l[ 
are disjoint. 

• Let u = E\eie 2 \ when h = I 2 , I'l and l^ share exactly one end¬ 
point, and I'l U I 2 is disjoint from li = I 2 . 

• Let V = E\eie 2 \ when li and I 2 share exactly one endpoint, l'^ 
and I 2 share exactly one endpoint, and Zi U I 2 and I'l U I 2 are 
disjoint. 

Then, q = s -\- 2{u -|- u) > 0. 

Arsuaga et al [I] use the above lemma to prove that the mean squared 
linking number of two uniform random polygons of length n is ^n‘^q where 
q is defined in Case (3) of the lemma. We now apply the above lemma in a 
similar way to obtain a formula for the mean squared linking number of two 
uniform random polygons where the number of vertices in the two polygons 
may differ. 
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Lemma 2.3. Let n, m > 3, and let the graph G he the disjoint union of an 
n-cycle L and an m-cycle L'. Then the mean squared linking number of a 
random linear embedding of G in the cube is ^nmq, where q is defined 
in Case (3) of Lem,ma \2.Si 

Proof. Let the edges of Lhe li,l 2 , ■ ■ ■, In and the edges of L' be l[, I 2 , ■ ■ ■, I'mi 
both cyclically ordered and oriented. Let Cij denote the signed crossing of 
li and I'j. Then, the linking number of L and L' is given by 

^ n m 

lk{L,L') = -J2Y.^^r 

i=l j=l 

Hence the expected value of the mean squared linking number is given 
by: 


E 


^ n m \ 

, i-ii-i / 
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Note that those cross terms which we know are 0 by Cases (1) and (2) of 
Lemma 12.21 have been omitted from the above expansion. 

Let 2s denote the probability that a pair of edges k and /'• cross. Then 
the first term in the above expansion is 


1 

4 


n m 


EE ^141 


1 

—nms. 
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Let u denote the expected value of the product of signed crossings of an 
edge li of L with consecutive edges /'■ and of L'. Then the second term 
in the above expansion is given by 

^ n m 

- + ^Ne*(i+i)]) = nmu. 

i=l j=l 

Let V denote the expected value of the product of signed crossings of con¬ 
secutive edges li and li±i of L with consecutive edges /'• and ll_^_l respectively 
of Lh Then the third term in the above expansion is given by 


^ n m 

-EE(-eN'( 

i=i j=i 
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Finally, let q = s + 2{u + v). Then we have 


E 


^ n m 
1=1 j=i 


2-1 




1 1 / / NX 1 

= —nms + nmu+nmv = —nm[s + 2[u+v)) = —nmq. 


□ 


Using the above lemma, we now prove the following theorem. 

Theorem 2.4. Let n > 6, and let Kn be a random linear embedding of the 
eomplete graph on n vertices in the cube . Then the mean sum of squared 
linking numbers for Kn is of the order of 6{n{n\)). 

Proof. Let k, I >3 such that k + l < n. k ^ I, then the number of disjoint 
pairs of cycles in Kn such that one cycle has k vertices and the other cycle 
has I vertices is given by 

fn\ /n — k\ {k — 1)! (/ — 1)! 

V/tyV ^ ) 2 2 ■ 

If /c = /, then this number is given by 

l/n\ l'n-k\ {k-l)\{l-iy. 

2VfeyV ^ ) 2 2 ’ 

By Lemma 12.31 we know that the mean squared linking number of a k- 
cycle and an /-cycle in Kn is where q is defined in Case (3) of Lemma 
ESI Thus, we obtain the mean sum of squared linking numbers over all 
disjoint pairs of cycles in Kn as 


!eeK:)("7')^^ 

fc=3 /=3 ^ ^ ^ ^ 



n—3 n—k 

EE 

k=3 1=3 


n — 
I 


k\ kl ll 
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fL 

16 


n—3 n—k 


EE 


nl 

(n-k-l)'.' 


Observe that the double sum 

n—3 n—k 


EE 


nl 

{n — k — l)\ 


counts the number of ways to obtain disjoint subsets of > 3 and I > 3 
ordered points from the set of n points. This same quantity can alternatively 
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be counted by choosing an ordered list ol i = k + I points out of n, then 
picking a number 3 < j < i — 3, so that the first j points are in the subset 
of k points, and the rest are in the subset of I points. Hence, we have the 
equality: 


n\ 


n—Sn—k 




If we only consider the i = n term in the sum, we obtain the following 
lower bound for the mean sum of squared linking numbers. 


> --—(z — 5) >- -In — 5 = —(n — 5)n!. 


16 ^ in — i)\ 
^=6 


For sufficiently large n, we have n —5 > S. Thus we have the lower bound 


i=6 ^ ' 


For an upper bound, we find that. 


2=6 ' 


n 


16 “ (n — i)\ 

2=6 ^ 

n 


^ (n-i)! 


2=6 

CXD 


< —n{n\) ^ 

16 ^ ml 

m=l 

9 ? n 

= —n(n!)e. 

16 


Putting these inequalities together, we see that the mean sum of squared 
linking numbers is of the order of 0(n(n!)). □ 


Another way to model entanglement is to consider the tangling of indi¬ 
vidual cycles rather than the linking between cycles. In particular, given 
a fixed oriented fc-cycle Jk in M^, we define the directional writhe Wr^(Jfc) 
projected in a direction perpendicular to a unit vector ^ G 5^ as the alge¬ 
braic sum of the signed crossings of J^. In order to avoid issues of sign, it is 
preferable to work instead with the direetional squared writhe, which is de¬ 
fined as Wr|(Jfc) = (Wr^(Jfc))^. Now if we average the directional squared 
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writhe Wr|(Jfc) over all possible direction vectors ^ G 5^, we obtain the 
mean squared writhe denoted by Wr^(Jfc). More formally, we define 

Wr2(Jfc) = ^ [ Wr|(Jfc)d^ 
dvr Js2 

Panagiotou et al m prove that the mean squared writhe of a random 
linear embedding of a /c-cycle is of the order of 0{k‘^). Rather than focusing 
on a single cycle, we are interested in obtaining a single value representing 
the complexity of the entanglement of all cycles C in a random linear em¬ 
bedding of Kn- Thus we define the mean sum of squared writhe of Kn as 
the expected value 


E[ ^ Wr2(C)] 

CCK„ 

over all random linear embeddings of Kn- 

We will make use of the following lemma from m which is similar to 
Lemma 12.21 

Lemma 2.5 f [T3]L Let denote edges in a random linear embed¬ 

ding of a graph with an orientation assigned to eaeh edge, let e* denote the 
signed erossing of h and l[, and let E[eie 2 \ denote the expeeted value of exe^- 
Also, let s, u, and v be defined in Case (3) of Lemma \2. 21 and letw = E[eie 2 \ 
when h, hj'i, I 2 o.re eonseeutive edges. Then q' = 2>s -\- 2(2u -|- n -|- rc) > 0. 

Theorem 2.6. Let n > 3, and let Kn be a random linear embedding of the 
eomplete graph on n vertices in the eube . Then the mean sum of squared 
writhe for Kn is of the order of 6{n{nl)). 

Proof. For some k < n, let Jk be a /c-cycle in Kn- It follows from Panagiotou 
et al m that the mean squared writhe satisfies 

Wr^(Jfc) = qk^ — {6q — q')k 

where q is defined in Case (3) of Lemma 12.21 Also, by Lemma 12.51 we know 
that q' > Q. Thus we have 


Wr2(Jfc) > qk^ -6qk. 


Hence, we have the following lower bound for the mean sum of squared 
writhe. 


n 

E[ Y, Wr2(J)] > 
JCKn k=3 



n\ 

k)\{2k) 


qn\ ^ /c — 6 

2 ^ in — k)\ 
k=3 ^ ' 
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Taking only the term when k = n, we see that 


qn\ ^ k — 6 ^ 

2 ^ in — k)\ ~ 
k=3 ^ ' 



6 ). 


For sufficiently large n, we have n — 6 > Thus we obtain the lower 
bound 


E[ Wt\J)] > (I) (n)n!. 

JCKn 

In order to get an upper bound, first observe that for any k-cycle Jk, by 
Panagiotou et al m we have 


Wr^(Jfc) = qk^ — (bg — q')k 

Thus, taken over all cycles J in we have the expected value 

n I 

- ( 6 , - 


Now, by Lemma 12.21 q > 0. Thus we obtain the following upper bound. 


- (6« - q')k) — 

Li (" 


ni 


k=3 


ni 


k)\{2k) ^'{n-h)\{2) 


n! qk + q' 

2 ^ (n — k)\ 
k=3 ^ ' 


I °° 1 

n\ 1 


, f . IL, ^ -I 

< {qn + qn}— 2^- 


< 


q + q' 


i=i 


(n)n!e. 


Putting these inequalities together we see, that the mean sum of squared 
writhe of a random linear embedding of Kn is of order of 0(n(n!)). 


□ 


We remark that the calculations above can be modified to show that the 
total number of links (resp. cycles) in a random linear embedding of iF„ is of 
the order of 0 (—), so that the mean average linking number (resp. writhe), 
where we average over all two component links (resp. cycles) in K^, is of 
the order of 9{n^). This agrees with the results of [T] and m, and shows 
that links (resp. cycles) of length 9{n) dominate the mean sum of linking 
number (resp. writhe) of the embedded graph. 
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3. Random linear embeddings of (n,p)-GRAPHS 

Theorem 3.1. For any knot J, the probability that a random linear embed¬ 
ding of an {n,p)-graph contains a cycle isotopic to J goes to 1 as n ^ oo. 

Proof. By Negami there is an integer R{J) such that every linear em¬ 
bedding of the complete graph Kr(j) in 1^^ contains J. Thus, given a set 
of vertices {vi,V 2 , ■ ■ ■ ,p_r(j)} in general position in the cube the linear 
embedding of Kr{j) defined by these vertices necessarily contains the knot 
J as a cycle. Furthermore, since no cycle in Kr(j) has length more than 
R{J), a random linear embedding of an (R(J),p)-graph has probability at 
least p^^R of containing the knot J. 

For any n > R{J), by partitioning the vertices into sets of R{J) vertices, 
we see that the probability that a random linear embedding of an {n,p)- 
graph contains J is at least 1 — (1 — . This value goes to 1 as 

n —>■ oo. □ 

Similar results can be obtained for other intrinsic properties of spatial 
embeddings of graphs. For example. 

Theorem 3.2. For an {n,p)-graph G, 

( 1 ) The probability that a random linear embedding of G contains a non¬ 
trivial link of two components is at least 1 — (1 — . In partic¬ 

ular, it goes to 1 as n ^ oo. 

(2) The probability that a random linear embedding of G contains a non¬ 
split link of three components is at least 1 — (1 — p®)!"/®!. In partic¬ 
ular, it goes to 1 as n —>■ oo. 

Proof. The proof is similar to that of Theorem 13.11 For part (1), we apply 
the result of Conway and Gordon [3] that every embedding of Kq contains a 
non-trivial link; and for part (2), we apply the result of Naimi and Pavelescu 
[9] that every linear embedding of Kg contains a non-split link of three 
components. □ 

Suppose that n > 6, p & (0) !)> and G is a random linear embedding of an 
(n,p)-graph. Then for any k, I > 3 such that k-\-l < n, the probability that 
G contains a pair of disjoint cycles where one is a /c-cycle and the other is 
an /-cycle is p^~^K We can now modify the proof of Theorem 12.41 to obtain 
the following. 

Theorem 3.3. For any n > 6 and p G (0,1), the mean sum of squared 
linking numbers of a random linear embedding of an (n,p)-graph is of the 
order of 6{p^n{n\)). 

Proof. We note that mean sum of squared linking numbers for a random 
linear embedding of an (n,p)-graph is 
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Taking the term i = n gives the lower bound: 

-^p”(n!)(n - 5) > ^p”n(n!). 
lb 62 

For the upper bound, we re-index over k = n — i, so that the sum becomes 

—5) (n — k — b) 

n\ > --^ = n\y — 

^ IT) — 7.V 


(n — i)! 

i =6 ^ ' k=0 


k\ 


n—6 

< n(n!K ^ 


/c=0 

CO _h 

p ^ 


kl 


< n(n!)p"' ^ 
k=0 


□ 


4. Random linear embeddings of Kq and iF3,3,i 


In this section, we apply the formulas in the proof of Theorem 12.41 to the 
graphs Kq and ^ 3 , 3,1 in order to find bounds on the probability of specific 
types of linking occurring. 

It follows from the formula 

16 ^ ^ in — k — ly. 

for the mean sum of squared linking numbers for iF„ in the proof of Theorem 
[M that the mean sum of squared linking numbers for Kq is given by 


Teo! 


45g. 


It has been shown independently by Hughes [7], Huh and Jeon [8], and 
Nikkuni [12] that every linear embedding of Kq contains either exactly one 
or three Hopf links, and all other links in the embedding are trivial. Hence, 
for a given linear embedding of Kq, the sum of squared linking numbers is 
either 1 or 3. 

This means that ibq = pi + 6pQ , where pi is the probability that a random 
linear embedding of Kq has exactly one Hopf link, and p 3 = 1 — pi is the 
probability that a random linear embedding of Kq has exactly three Hopf 
links. This implies that 

3 -45u 


Pi = 


2 
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We estimated the value of q to be 0.033867 ±0.000013 (see the numerical 
computation described in Appendix (21). This value is consistent with the 
value oi q = 0.0338 ± 0.024 obtained by [T], which is verified in [13]. Thus 
we find that pi = 0.7380 ± 0.0003. 

For the complete tripartite graph Naimi and Pavelescu [TO] show 

that every linear embedding contains either 1,2,3,4, or 5 non-trivial links. 
Furthermore, they show that if the number of non-trivial links is odd, all 
such links are Hopf links; whereas if the number of non-trivial links is even, 
then one link is a (2,4)-torus link and the rest are Hopf links. Since a Hopf 
link has linking number ±1 and a (2,4)-torus link has linking number ±2, 
it follows that the sum of squared linking numbers for any linear embedding 
of is either 1, 3, 5 or 7. 

Now every pair of disjoint cycles in it's, 3,1 consists of one 3-cycle and one 
4-cycle. By Lemma [2.3l the mean squared linking number of a random linear 
embedding of a disjoint union of a 3-cycle and a 4-cycle is |(3)(4). Since 
there are nine pairs of disjoint cycles in it follows that the expected 

value of the sum of squared linking numbers of a linear embedding of 
is 


|(9)(3)(4) =54g. 

For each k, we let pk be the probability that there are k non-trivial links 
in the embedding. Then, this expected value is equal to 

Ipi ± 5p2 ± 3p3 ± 7p4 ± 5p5 > Pi ± 3(1 - pi). 

Hence, it follows that the probability that there is precisely one non-trivial 
link in a random linear embedding of ^^ 3 ^ 3 ^! is given by 

3 - 54g 

Pi >- = 0.5856 ± 0.0004. 

5. Experimental data 

In this section, we describe some experimental results we obtained for 
links in random linear embeddings of graphs. 

The data was generated using a Python program, taking coordinates of 
the n vertices to be uniformly distributed in (0,1). An edge between two 
vertices is taken with probability p, and then the number of links with 
each linking number are tallied. To give a more accurate picture of the 
distribution of linking numbers and the average sum of linking numbers, we 
took multiple samples for each (n,p). 

We first investigated the mean sum of squared linking number for p = 1, 
p = 0.5, and p = 0.25, comparing experimental data with the expected 
number from the formula given in the proof of Theorem 13.31 For p = 1, we 
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Figure 1. Experimental vs. expected mean sum of squared 
linking number for p = 1. Experimental data used 1000 
samples for 6 < n < 11, and 100 samples for n = 12. 
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Figure 2. Experimental vs. expected mean sum of squared 
linking number for p = 0.5. Experimental data used 1000 
samples for 6 < n < 14, and 100 samples for n = 15. 
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Figure 3. Experimental vs. expected mean sum of squared 
linking number for p = 0.25. Experimental data used 200000 
samples for 6 < n < 7, 50000 samples for 8 < n < 12, 5000 
samples for 13 < n < 15, 500 samples for n = 16, and 100 
samples for n = 17 . 


took 1000 samples for 6 < n < 11, and 100 samples for n = 12 (see Eigure 
[1]) For p = 0.5, we took 1000 samples for 6 < n < 14 and 100 samples for 
n = 15 (see Figure [2]). For p = 0.25, we took 200,000 samples for 6 < n < 7, 
50,000 samples for 8 < n < 12, 5000 samples for 13 < n < 15, 500 samples 
for n = 16, and 100 samples for n = 17 (see Figure [3]). 

The experimental data follows the expected super-factorial growth. The 
deviation for the mean sum of squared linking number is within approxi¬ 
mately 10% of the expected value (and most data points are within 5%), 
and the discrepancy for large numbers of vertices and p < 1 is due to the 
small number of samples taken due to computational constraints. In addi¬ 
tion, from the n = 6, p = 1 case, we can determine that of the 1000 random 
linear embeddings sampled, 729 had exactly one Hopf link, giving a 99% 
confidence interval for the probability that a random linear embedding of 
Kq has one Hopf link of 0.729 ± 0.036, which agrees with the theoretical 
computation from Section [4] and numerical value of q. 

In addition, we computed the mean average squared linking number and 
the mean average absolute linking number, where the average is taken over 
all links in the graph, and then the mean is taken over all samples of a given 
size. The experimental data for the mean average squared linking number 
follows a quadratic growth, as expected. From the samples we computed, 
it appeared that the average absolute linking number was also quadratic. 
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Mean average squared and absolute linking number 



Figure 4. Mean average squared and absolute linking num¬ 
ber for p = 1. 

Mean average squared and absolute linking number 



Figure 5. Mean average squared and absolute linking num¬ 
ber for p = 0.5. 


which differs from the linear growth rate for absolute linking number of 
random polygons studied in [I] and m (see Figures miSl and [6]). However, 
because of the small number of data points that we could compute, this is 
inconclusive. 
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Mean average squared and absolute linking number 



Figure 6. Mean average squared and absolute linking num¬ 
ber for p = 0.25. 



Figure 7. Proportion of links in random linear embeddings 
of Kn with linking number 0, 1, and 2. 

We also investigated the distributions of links with a given linking number 
in random (n,p) graphs (see Figures [7l [8l and [9]). Omitted from the figures 
are links with linking number greater than 2, which were detected in fewer 
than 1% of the links in the samples that we generated. We expect that for a 
hxed linking number k, the proportion of links with that linking number will 
increase, peak, then decrease as n —oo. However, due to computational 
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Figure 8 . Proportion of links in random linear embeddings 
of (n,p) graphs with linking number 0, 1, and 2, when p = 
0.5. 



Figure 9. Proportion of links in random linear embeddings 
of (n,p) graphs with linking number 0, 1, and 2, when p = 


0.25. 


constraints, we were not able to compute samples out to large enough n to 
see this behavior, even with linking number 1. 
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Appendix A. Computing q 

In order to obtain better estimates of the linking probabilities for Kq 
and we numerically estimated the value of q. Consider two triangles 

described by the consistently oriented edges h and /i, ^ 2 ) ^ 3 ) and let eij 
denote the signed crossing number between li and I'j. Then it follows from 
the proof of Lemma 12.21 in [1] that 


E 



18q. 


But the quantity Ylij=i precisely twice the linking number of the 
two triangles, which for a linear embedding is either 0 or ±1. Hence, ^ is 
the probability that a random linear embedding of two disjoint triangles is 
linked. 

We wrote a Python program to generate random linear embeddings of two 
triangles by taking six random points in C^, described as three (pseudo)random 
coordinates in (0,1), and computing the associated linking number. Out of 
the one billion samples generated, 152,402,780 were linked, giving a 99% 
confidence interval for q of 0.033867 ± 0.000013. 

All code is available at http: //math. berkeley. edu/~kozai/random_graphs/, 
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